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ABSTRACT 

The combination of mid- infrared and recent submm/mm measurements allows 
us to set up the first comprehensive spectral energy distribution (SED) of the 
circumstellar material around a young Brown Dwarf. Simple arguments suggest 
that the dust is distributed in the form of a disk. We compare basic models to 
explore the disk parameters. The modeling shows that a flat disk geometry fits 
well the observations. A flared disk explains the SED only if it has a puffed- 
up inner rim and an inner gap much larger than the dust sublimation radius. 
Similarities and differences with disks around T Tauri stars are discussed. 

Subject headings: circumstellar matter, stars: low-mass, brown dwarfs, stars: 
individual (CFHT-BD-Tau 4) 



1. Introduction 

Hundreds of brown dwarfs (BDs) have been discovered in the last few years, both within 
star-forming regions (e.g. Moraux et al. 2003; Persi et al. 2000) and among field stars (e.g. 
Phan-Bao et al. 2003; Kirkpatrick et al. 1997). Still very little is known about their formation 
processes. 

First indications that BDs may form as ordinary stars, i.e. via accretion from a cir- 
cumstellar disk, come from observations of Ha ( Jayawardhana et al. 2003) and near-infrared 
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excess emission (Muench et al. 2001; Liu et al. 2003). However, in the case of BDs only 
excess emission at longer wavelengths provides the indisputable evidence for the presence of 
circumstellar material (Natta & Testi 2001). Thermal emission of warm dust around BDs 
has been observed by a number of groups (Comeron et al. 2000; Persi et al. 2000; Testi et 
al. 2001; Apai et al. 2002). Very recently, our search for submm/mm emission from young 
BDs provided the first detection of cold dust around two of them (Klein et al. 2003). For 
the first time, this allowed to estimate the mass of the circumstellar material and to relate 
it to the mass and luminosity of the central object. 

These results open up new interesting questions: Is the circumstellar material around 
BDs confined in a disk? Are these disks similar to T Tauri disks or different ones? What 
is the effect of the ultra-low luminosity on the disk properties? Tentative modeling of disks 
has been performed by Natta & Testi (2001) based on two ISOCAM mid-infrared (MIR) 
measurements with the general assumption of flared disk geometry. The first ground-based 
observations of one of their objects showed that the flared disk model cannot be applied 
blindly: at least in one case a flat configuration can better explain the observations (Apai 
et al. 2002). The same flat geometry has been later used by Natta et al. (2002) to interpret 
the ISOCAM measurements of other BDs. 

A thorough investigation of the disk properties requires measurements of different disk 
zones, i.e. in different wavelength domains. In this Letter we compile the first complete 
spectral energy distribution (SED) of a disk around a BD. The goal of this paper is to 
test the basic models and determine which model parameters can be constrained from the 
combination of short- and long-wavelengths observations. Our study focuses on the object 
CFHT-BD-Tau 4 (2MASSI J0439474+260140) which was found and confirmed as a young 
BD by Martin et al. (2001). Its very strong Ha emission (Martin et al. 2001) coupled with 
X-ray activity (Mokler & Stelzer 2002) and further spectral signatures (Jayawardhana et al. 
2003) prove that this BD is intensively accreting matter. 

In the following, we first present the SED of CFHT-BD-Tau 4, we prove that the cir- 
cumstellar material cannot be in a spherical shell and we argue towards a disk configuration. 
We model the SED applying different disk geometries. Finally, we discuss similarities and 
differences with accretion disks around T Tauri stars. 



2. The Spectral Energy Distribution 

In this section we compile the observations available on CFHT-BD-Tau 4. Optical and 
near-infrared data, together with an estimate of the line-of-sight extinction (A v =3 mag), 
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have been published by Martin et al. (2001). An additional L-band flux is given in Liu et al. 
(2003). We retrieved archival ISOCAM images at 6.7 and 14.3 yum and we measured fluxes 
of 28.4 and 38.5 mJy by using an aperture of 24". We accept a photometric accuracy of 20% 
as given in the ISOCAM Calibration Accuracies Document 1 . Our group observed CFHT- 
BD-Tau 4 in the mm and submm continuum at the IRAM-30m and the JCMT telescopes. 
The data reduction and analysis is presented in Klein et al. (2003). We provide the summary 
of the above fluxes and errors in Table 1. 



3. Input Parameters 

The main input parameters for the modeling are derived directly from the observations. 
CFHT-BD-Tau 4 is expected to be a member of the Taurus star-forming region located at 
a distance of 140 pc. We assume an effective temperature of 2900 ± 100 K for an M7±0.5 
spectral type (Martin et al. 2001; Briceno et al. 2002) and a mass of 70 M Jup . In agreement 
with Martin et al. (2001), we calculate a total luminosity of 0.1 L which is brighter than 
the luminosity of the youngest (1 Myr) isochrone from Allard et al. (1999). 

The wavelength dependence of the dust opacity (k\) and the mass of the circumstellar 
dust can be calculated from the expected optically thin submm and mm fluxes. In the 
optically thin regime the flux depends on the wavelengths as oc A~( 2+/3 \ where (3 is the 
wavelength dependence of k x . A value of (3 — 2 has been measured for interstellar dust 
particles (Draine & Lee 1984), while (3 — 1 is generally used for circumstellar disks where 
grain growth is expected to take place (Beckwith et al. 2000). In our case the slope of the 
SED is —3.8 ± 0.9, with the error computed from the uncertainty in the flux measurements. 
This implies that (3 is 1.8 ±0.9. Due to the large error on the measured (3, we cannot fix the 
dust opacity. To compare our results with disks around T Tauri stars, we assume the same 
dust opacity as suggested by Beckwith et al. (1990) which amounts to 2 cm 2 /g at 1.3 mm. 
For a dust temperature between 10 and 20 K and a gas-to-dust ratio of 100, we calculate a 
total circumstellar mass of 1.5-0.3 M Jup . We note that adopting the 10 times lower opacity 
of interstellar grains would lead to 10 times larger circumstellar masses. 



1 http://www.iso.vilspa.esa.es/manuals/CAM/accuracies/CAM J accur.pdf 
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4. Modeling 

A simple spherical shell model is first investigated by using the radiative transfer (RT) 
code of Manske et al. (1998). The mass of the circumstellar material is set to 4 Mj up in 
order to reproduce the submm and mm fluxes. Different power-law density distributions, 
inner/outer radii and dust properties were tested (see Table 2). The ID simulations always 
predict a visual extinction at least an order of magnitude higher than measured by Martin 
et al. (2001) and Itoh et al. (2002). Thus, the simple shell configuration has to be excluded 
and a 2D distribution has to be considered. 

The most plausible 2D circumstellar dust distribution is believed to be in the form 
of a disk (Terebey et al. 1993). In the following, we test various versions of the most 
widely accepted disk configurations, i.e. the flat (Adams et al. 1987) and the flared disk 
(Kenyon & Hartmann 1987). The opening angle of the flat disk is constant, while that of the 
flared disk increases with radius. An optically thin superheated layer is expected for flared 
configurations. An addition to some of the flared disks is the self-consistent treatment of the 
inner rim that, due to its higher temperature, puffs-up with respect to the flaring disk behind 
it (Dullemond et al. 2001). Similar changes in the flared disk geometry have been proposed 
by Chiang & Goldreich (1997) as an effect of stellar accretion and viscous dissipation. For 
all the models, we adopt a dust opacity which follows the silicate mass absorption coefficient 
by Draine & Lee (1984) for short wavelengths and the opacity law by Beckwith et al. (1990) 
for long wavelengths. The transition wavelength depends on the grain size: for small grains 
of 0.1 ^m radius the two opacity laws meet at 100 ^m, while for 10 fxm grain they meet at 
200 urn. 

We consider two flat disk models (m4 and m5 in Table 2): one with a single disk 
temperature distribution (Beckwith et al. 1990) and one with an additional hot surface 
layer. The main difference between the two approximations is at 10 /txm, where the two-layer 
model predicts an emission feature, which is absent in the optically thick one-layer flat disk. 
For both configurations, plausible physical parameters provide a good fit. Because of the 
lack of measurements at the silicate feature, we cannot discriminate between the two cases. 

For the flared disk configurations (Dullemond et al. 2001), a too high MIR flux is pre- 
dicted by the models when the disk inner radius is determined by the dust sublimation 
temperature. To reduce the MIR emission three possibilities are considered: high disk in- 
clination, a large inner gap and the shadow of a puffed-up inner rim. Considering the low 
line-of-sight extinction measured towards CFHT-BD-Tau 4, almost edge-on disks (inclina- 
tions larger than 60°) are not plausible. Even when a large inner gap of about 0.07 AU is 
introduced in the classical flared disk model (m7), the MIR measurements at 6.7 and 14.3 fim 
cannot be simultaneously fitted. On the contrary, the flared disk with a puffed-up inner rim 
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and large inner gap reproduces well the SED of CFHT-BD-Tau 4 (m9 in Table 2 and solid 
line in Fig. 1). 

To demonstrate the effect of grain growth, we also tested an opacity law for silicate 
grains with 10 fim radius (Draine & Lee 1984). Such grains lack the emission feature at 
A = 10 /im. The one-layer flat disk model is not affected by the different dust opacities. For 
the two-layer flat configuration, the silicate feature at 10 /im vanishes, but the SED can still 
be well fitted with parameters similar to those provided in Table 2. Because larger grains 
have lower opacities in the MIR regime, the flared disk model with the puffed-up rim can 
fit the observations with a slightly smaller inner gap (0.04 AU). Even with the addition of 
larger grains, the classical flared disk model does not explain the observed SED. 

5. Discussion 

The combination of short and long wavelengths observations allow us, for the first time, 
to investigate the properties of the circumstellar dust around a BD. The SED shows an 
overall shape similar to that of T Tauri stars. Adopting the opacity by Beckwith et al. 
(1990), we calculate a circumstellar mass of 0.3-1.5 Mj up , which is a few percent of the 
BD's mass. A similar mass ratio is found for T Tauri stars and their disks. Thus, the disk 
is unlikely to be truncated by any previous ejection event (Reipurth & Clarke 2001). A 
discussion is presented in Klein et al. (2003). 

Since spherically distributed dust would result in too high extinction, we argue towards 
the presence of a circumstellar disk. Our disk models suggest that either the disk is flat, or 
it is flared with a puffed-up inner rim and a large inner gap. A flat geometry was previously 
suggested by Apai et al. (2002) to explain the SED of a disk around a BD candidate. The 
present results support a flat disk as well, thus suggesting that this might be general for 
disks around very low-luminosity objects. A flat dust distribution could be a hint towards 
large grain sizes since large grains are expected to settle in the disk midplane (see Kenyon & 
Hartmann (1987) and Chiang & Goldreich (1997) for a discussion of the physical conditions 
under which flared disks can exist). Future observations at 10 /im could further strengthen 
this picture if no silicate feature is found. 

Should the disk be flared with a puffed-up inner rim, an inner gap 5 times larger than 
the dust sublimation radius is required to account for the observed SED. The presence of 
magnetic fields is not sufficient to explain inner gaps as large as those found in the modeling 
(about 9 R ), see for example Konigl (1991). The material close to the BD could be cleared 
out by a massive planet, as suggested in the case of TW Hya by Calvet et al. (2002). 
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However, the most realistic explanation is the presence of a close companion with a spectral 
type similar to the primary. This would also explain the over-luminosity of this BD (Martin 
et al. 2001). Unfortunately, the estimated separation of the possible companion would be 
< 1 mas, unresolvable with any current imaging technique. 

To fully characterize BD disks, a detailed exploration of their SED is essential. Es- 
pecially important are the silicate emission features and the wavelength range between 60- 
180 //m. The upcoming SIRTF and SOFIA missions will be capable of gathering a statisti- 
cally relevant number of complete SEDs, as well as providing a reliable disk frequency. 

6. Conclusions 

We present the first comprehensive models of the circumstellar dust surrounding a young 
BD. Our results are: 

• the overall shape of the SED is very similar to those of T Tauri disks 

• a spherical dust distribution contradicts the observed extinction 

• a simple flat disk fits well the observations 

• the flat geometry hints towards large grains 

• the classical flared disk does not reproduce the SED of CFHT-BD-Tau 4 

• a flared disk with a puffed-up inner rim can only explain the observed SED if it has 
an inner gap much larger than the dust sublimation radius 
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Table 1. Observed fluxes for CFHT-BDTau 4 



A 




Error Reference 


H 


[mJy] 


[mJy] 




0.658 


0.062 


0.002 


1 


0.822 


1.23 


0.03 


1 


1.25 


22.9 


1.0 


1 


1.65 


38.7 


1.0 


1 


2.17 


46.5 


1.4 


1 


3.75 


52.5 


1.5 


2 


6.7 


28.4 


5.7 


3 


14.3 


38.5 


7.7 


3 


850 


10.8 


3.7 a 


4 


1300 


2.14 


0.7 a 


4 


a The 


large 


uncertainty 


in the 



submm and mm measurements is due 
to the uncertain flux of the calibra- 
tors. The estimated statistical errors 
are smaller (1.8 and 0.56 mJy, respec- 
tively). 

References. - (1) Martin et al. 
(2001);(2) Liu et al. (2003);(3) This 
paper; (4) Klein et al. (2003) 
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Table 2. Summary of the Model Parameters and Results 



ID 


Type a 


Geom. 


Par am. b 


Rin 

[AU] 


Fits? c 


M disk 
[M Jup ] 


Comments d 


Refs. 


ml 


RT 


Shell 







No 


4 


A v ~60 mag 


(1) 


m2 


RT 


Shell 


-1 




No 


4 


A v ~300 mag 


(1) 


m3 


RT 


Shell 


-2 




No 


4 


A v ~1000 mag 


(1) 


m4 


AA 


Flat 


one-layer 


0.02 


Yes 


1.4 


i=20-45°(35°) 


(2) 


m5 


AA 


Flat 


two-layer 


0.01 


Yes 


4 


2 =5-30°(10°) 




m6 


sec 


Flared 


No Rim 


0.03 


(Yes) 


0.7 


i=80-82° 


(3) 


m7 


sec 


Flared 


No Rim 


0.07 


No 


0.7 




(3) 


m8 


sec 


Flared 


Rim 


0.03 


No 


0.3 




(3) 


m9 


sec 


Flared 


Rim 


0.05 


Yes 


0.3 


?=8-20°(12°) 


(3) 



Note. - Summary of the models and results. Given are the models' identifier, 
the geometry (shell or disk type), the characteristic parameters, inner disk radius, 
possibility of fitting the observations, total disk mass, comments and references 

a RT - Radiative Transfer Code; AA - Analytical Approximation; SCC - Self- 
Consistent Calculation 

b For ID models we provide the exponent of the power-law density distribution 
C A parenthesis stands for models which are physically not plausible 

d A v stands for visual extinction, % is the disk inclination (0°for face-on) with the 
best disk inclination in parenthesis 



References. — (1) Manske et al. (1998); (2) Beckwith et al. (1990); (3) Dullemond et 
al. (2001) 
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SED of CFHT-BD-Tau 4 




0.1 1.0 10.0 100.0 1000.0 10000.0 

X [am] 

Fig. 1. — Comparison of model fits with observations. The measurements are de-reddened 
with a visual extinction of 3 mag (Martin et al. 2001) and are plotted with asterisks. The 
photometric errors are taken from Table 1. The models that best fit the observations are the 
flat disk (dashed line, model m4 in Table 2) and the flared disk with an inner puffed-up rim 
(solid line, model m9 in Table 2). The atmosphere of a 1 Myr old BD with 3000 K effective 
temperature is overplotted in dots (Allard et al. 1999). CFHT-BD-Tau 4 is more luminous 
than the youngest isochrone from Allard et al. (1999). 
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